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Abstract 

We apply the statistical measure of complexity introduced by Lopez- 



o 

i-^ • Ruiz, Mancini and Calbet [1] to neutron stars structure. Neutron 

c/3 , stars is a classical example where the gravitational field and quantum 

^ ' behavior are combined and produce a macroscopic dense object. Ac- 

tually, we continue the recent application of Sahudo and Pacheco [2] 
pi . to white dwarfs structure. We concentrate our study on the connec- 

Tlj" I tion between complexity and neutron star properties, like maximum 

' mass and the corresponding radius, applying a specific set of realistic 

. equation of states. Moreover, the effect of the strength of the gravi- 

jy-^ I tational field on the neutron star structure and consequently on the 

' complexity measure is also investigated. It is seen that neutron stars, 

■ consistent with astronomical observations so far, are ordered systems 

(low complexity), which cannot grow in complexity as their mass in- 
creases. This is a result of the interplay of gravity, the short-range 
^ ' nuclear force and the very short-range weak interaction. 

H ; 

5r ! Keywords: Shannon Entropy; Complexity; Self-Organization; Equation of 

state; Neutron stars. 



1 Introduction 

Information theory, founded by Shannon to provide a theoretical framework 
in communications [3], has been further employed as a useful tool to charac- 
terize physical systems during the next decades [4, 5, 6, 7, 8]. Since then, a 
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series of studies concerning the application of information theory to various 
physical systems has appeared in the literature, while their number exhibits a 
remarkable rise over the past decade. The dependence of various information- 
theoretic measures on some parameters of the physical systems has been 
studied, the presence of correlations has been quantified, connections with 
experimental data have been detected and universal properties have been 
extracted. Very recently, these investigations have been extended to include 
statistical complexity measures, in order to examine self-organizing charac- 
teristics of physical systems, patterns and correlations. Although a complete 
and universal definition for complexity is missing, the current framework 
provides interesting and satisfactory results expected from intuition. So far, 
various complexity measures are used taking into account conditions and con- 
straints imposed by the physical system under consideration [9, 10, 11, 12]. 
Such information and complexity studies (focusing so far on the two sta- 
tistical complexity measures SDL and LMC [1, 13, 14]) have been apphed 
to various quantum many-body systems i.e. nuclei, atoms, atomic clusters, 
bosons and molecules [15]- [33]. Recently Saniido and Pacheco [2] were the 
first to extend those studies to an astronomical object i.e. a white dwarf. 
Specifically, the Shannon information entropy S and the statistical complex- 
ity C have been calculated in the two kinematics extremes (non-relativistic 
and relativistic cases) of the electron-degenerate matter of a dwarf. 

In the present work we study the information content of another astro- 
nomical object, a neutron star. With a mass of 1.4 up to 3 solar masses, a 
radius of ~ 10 Km and an average density of 10^^ g/cm^, the neutron star 
is one of the possible endpoints of stellar evolution. Further gravitational 
collapse is counterbalanced by repulsive forces originating from Pauli's ex- 
clusion principle, if the mass of the compressed stellar core is less than the 
Oppenheimer-Volkoff limit of about 3 solar masses. 

Neutron stars are systems with several similarities with atomic systems, 
but there are also fundamental differences. In an atomic system self-organization 
is reached through the competition between the Coulomb interaction and 
the Pauli principle. In fact, the long-range electromagnetic interaction is the 
main interaction among the particles of the system. In addition, atoms are 
microscopic systems with a typical dimension of a few Angstroms (10~^°m). 
In contrast, a neutron star is a macroscopic system with typical dimension of 
10^ m, much more complicated than the atoms, in the sense that it is orga- 
nized under the competition of mainly the following three interactions. The 
long-range force of gravity, whose pressure tends to compress the mass of 
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the star. The short-range nuclear interaction, which through the degeneracy 
pressure of the nucleons tends to extend the outer mass of the star. Finally, 
the very short-range weak interaction, which is a kind of regulator of the par- 
ticle fraction and thus affects indirectly the properties of the neutron star. 
The above forces coexist in harmony in the interior of a neutron star. The 
main features of the structure of a neutron star (mass, pressure and radius) 
are described by the Tolman-Oppenheimer-Volkoff equations [34], while they 
also depend strongly on the applied nuclear equation of state. 

In this Letter we present a study of the information properties of a neutron 
star and explore how they are connected with the characteristic properties of 
the structure of the system, i.e. its mass M and radius R. Furthermore, we 
investigate the dependence of information and complexity measures on the 
nuclear forces, through the asymmetry energy parameter c, and the gravita- 
tional constant G. Also we comment on the effect of information measures 
on the stability of a neutron star, based on the fact that stability regions are 
characterized by the inequality dM/dR < 0. 

Here, we consider that the temperature of a star is T = 0, in the con- 
text that the Fermi energy is much greater than kT. However, it should 
be of interest to extend our study and try to connect the thermodynamic 
properties of a hot neutron start with the information content of the sys- 
tem. Furthermore, it is important to examine the information properties of 
other astronomical objects e.g. stars consisting of fermions or bosons, with 
arbitrary masses and interaction strengths. Such a work is in progress. 

The outline of this Letter is the following: In Section 2, we define the 
information and complexity measures employed here, together with a model 
of neutron stars. In Section 3, we present our results and a discussion, while 
Section 4 contains a summary. 



2.1 Theoretical information measures 

The Shannon information entropy S [3] for a continuous probability distribu- 
tion p(r), denoting a measure of the amount of uncertainty associated with 
a probability distribution, is defined as 



2 The model 
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while the disequihbrium D, being a quadratic distance from equiprobabihty, 
is 



with dimension of inverse volume. 

For a continuous probability distribution the disequilibrium is indeed the 
same measure as the information energy defined by Onicescu [35]. 

For discrete probability distributions {pi} = {pi,p2,--- ,Pn}, the infor- 
mation entropy S = —Yl!i=iPi ^^P* minimum (^min = 0) for the dis- 
tribution of a completely regular system (absolutely localized), where one 
of the Pi's equals unity, while all the others vanish. The maximum value 
(>S'max = InA?") is attained for the equiprobable distribution (completely de- 
localized), where pi = l/N,i = l,...,N. On the other hand, the disequilib- 
rium D = ^^i(l/A^ ~ PiY 1 is maximum, -Dmax = 1 — 1/-^ ~^ ^ (f*^^ large 
N) for a completely regular system, while it is minimum, Dmin = for an 
equiprobable distribution. 

In the continuous case, an equiprobable probability distribution can be 
defined as a rectangular function, while a completely regular system corre- 
sponds to a 5-like probability distribution function, where the width of the 
distribution becomes very narrow and its peak extremely high. 

In order to study the statistical complexity defined by Lopez-Ruiz, Mancini 
and Calbet (LMC) [1] , we use a slightly modified definition introduced in [36] 



is the information content of the system, while the exponential functional 
preserves the positivity of C. 

The aforementioned definitions of information entropy and disequilibrium 
in the case of neutron stars are modified as follows: 




(2) 



C^H-D, 



(3) 



where 




(4) 





and 




(6) 
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where Bq = 8.9 x 10"^ Km~^ is a proper constant satisfying the condition 
that both information entropy S and disequilibrium should be dimensionlcss 
quantities, while e(r) is the dimensionlcss energy density of the system. It 
is equivalent to the density mass p(r), obtained by solving the structure 
equations characterizing the system. 

2.2 Neutron star structure equations 

In order to calculate the gross properties of a neutron star, we assume that 
the star has a spherically symmetric distribution of mass in hydrostatic equi- 
librium and is extremely cold (T = 0). Effects of rotations and magnetic 
fields are neglected and the equilibrium configurations are obtained by solv- 
ing the Tolman-Oppenheimer-Volkoff (TOV) equations [34] 

GM{r)p{r) f P(r) \ f Anr'P{r) \ f _ 2GM(r) V' 
r2 V ^ cV(r) J \ c^M{r) )\ c^r ) ' 

47rrV r = 7 

where P(r) and M(r) are the pressure and the mass functions of the star 
respectively. 

To solve the set of equations (7) for P(r) and M(r). one can integrate 
outwards from the origin (r = 0) to the point r = R, where the pressure 
becomes zero. This point defines R as the radius of the star. To do this, one 
needs an initial value of the pressure at r = 0, called Pc = P{r = 0). The 
radius R and the total mass of the star, M = M{R), depend on the value of 
Pc- To be able to perform the integration, one also needs to know the energy 
density e(r) (or the density mass p(r)) in terms of the pressure P(r). This 
relationship is the equation of state (EOS) for neutron star matter and here, 
has been calculated applying a phenomenological nuclear model. 

We can modify equations (7), so that they become suitable for program- 



dP{r) 
dr 

dM{r) 
dr 
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ming, in the following form: 

dP(r) e(r)M(r) f P(r)\ ( Pir 

- -1.474 ^ ^ , ^ M 1 + 1 + 11.2 X 10~V^ ^ 



dr \ e(r) / \ M(r) 



X 1-2.948 



dM(r) 
dr 



11.2 X 10-^ e(r). (8) 



In Eqs. (8), the quantities P(r), e(r) and M(r) are dimensionless. The radius 
r is measured in Km. More specifically: 

M(r) = M(r) M©, e(r) = e(r) eo, P(r) = P(r) eo, eo = 1 MeV fm^l 

(9) 

It is obvious from Eqs. (8) that 

M{R) = 11.2 X 10"^ e(r)dr = bo J e(r) dr. (10) 

2.3 Nuclear equation of state 

In general, the energy per baryon of neutron-rich matter may be written to 
a very good approximation as 

^^^ + (l-2.f^.,„W, (11) 

where n is the baryon density (n = n„ + Up) and x is the proton fraction 
{x — np/n). The symmetry energy Egymin) can be expressed in terms of the 
difference of the energy per baryon between neutron {x = 0) and symmetrical 
{x = 1/2) matter. Here, we consider a schematic equation for symmetric 
nuclear matter energy (energy per baryon E/A or equivalently the energy 
density per nuclear density e/n), given by the expression [37] 

E(n, 1/2) ^e^^ ^ 3^, ^ ^ u^n/uo, (12) 

An 5 

where Ep = (3/5)(/lA:^)^/(2miv) is the mean kinetic energy per baryon in 
equilibrium state and no is the saturation density. 
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The density dependent potential V(u) of the symmetric nuclear matter 
is parameterized, based on the previous work of Prakash et. al. [37, 38], as 
follows 

where pp is the Fermi momentum, related to by pp — p%u^^^. The 
parameters Ai and A2 parameterize the finite-range forces between nuclcons. 
The values employed here are Ai = 1.5p^ and A2 = 3pp. The parameters A, 
B, B', a, Ci and C2 are determined using the constraints provided by the 
empirical properties of symmetric nuclear matter at the saturation density 
no- Then, the values of the above parameters are determined in order that 

E{n = no) /A - uinc^ = -16 MeV, Uq = 0.16 fm-^ Kq = 240 MeV. 

In general, the parameter values for three possible values of the compression 

modulus Kq (^Kq = 9nQ^^^^5^|„(, j arc displayed in Table 1, in [37]. 

To a very good approximation, the nuclear symmetry energy £^sym can be 
parameterized as follows [39] 

Esym(M) ^13ix'/=^ + 17F(m), (14) 

where the first term of the right-hand side part of Eq. (14) is the contribution 
of the kinetic energy and the second term comes from the interaction energy. 
For the function F{u) , that parametrizes the interaction part of the symmetry 
energy, we apply the following form 

F{u) = u", (15) 

where the parameter c (hereafter called potential parameter) varies between 
0.4 < c < 1.5 leading to reasonable values for the symmetry energy. In order 
to construct the nuclear equation of state, the expression of the pressure is 
needed. In general, the pressure, at temperature T = 0, is given by the 
relation 

„ , dit/n) de , 
P = ^' ' = n — - e. 16 
dn dn 

Employing equations (11), (12) and (16), we find the contribution of the 
baryon to the total pressure: 



2 po ^ ,,5/3 , dyju) 

5 du 



\2 „,2 dEsyTD.{u) 



+ no{l-2xYu' —^y^^-\ (17) 
du 
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The leptons (electrons and muons) , originating from the condition of the 
beta stable matter, contribute also to the total energy and total pressure [39]. 
To be more precise, the electrons and the muons, which are the ingredients 
of the neutron star, are considered as non-interacting Fermi gases. In that 
case their contribution to the total energy and pressure is 

= ^ + + - ' ^'^^ 

^ [^2"' " ^"^^^ + + 3sinh-^(.)] , (19) 

where z = kp/mic. Now the total energy and total pressure of charge neutral 
and chemically equilibrium nuclear matter are 

etot = Cft + XI 
Ptot^Pb+ Yl (21) 

l=e~,ix~ 

Prom equations (20) and (21) we can construct the equation of state in the 
form e = e(P). In order to calculate the global properties of the neutron star, 
i.e. the radius and mass, we solved numerically the TOV equations (7) with 
the given equations of state constructed employing the present model. For 
very low densities {n < 0.08 fm~^) we use the equation of state according 
to Feynman, Metropohs and Teller [40] and also from Baym, Bethe and 
Sutherland [41]. 



3 Results and Discussion 

The starting point of our study is the solution of Eq. (8) for three different 
equations of /3-stable nuclear matter. More precisely, we employ three values 
of the parameter c, which characterizes the density dependence of the nuclear 
symmetry energy, i.e. c = 0.7 (soft equation of state), c = 1.0, and c = 1.5 
(stiff equation of state). In Fig. 1(a), we plot the nuclear symmetry energy 
Esyia, in Fig. 1(b) the corresponding equations of state and in Fig. 1(c) the 
mass-radius diagrams for each of the three cases. 

Actually every pair {R, M) in a mass-radius diagram is the outcome of the 
structure equations (Eqs. 8) for an arbitrary chosen initial value of the pres- 
sure Pc in the center of the star. Thus, varying the value of Pc in a reasonable 
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range, we can have a picture of the behavior of those substantial structure 
characteristics. We have to note here that the region where dM/dR < 
corresponds to a stable neutron star, while dM/dR > to an unstable one. 
The presence of the unstable region (M < Mmax and dM/dR > 0) seems to 
lead to double valued functions S{M) and C{M), for values of M close to 
-^max (insets in Fig. 2(a) and Fig. 2(d)). However, the study of that region 
is beyond the scope of the present work. Another important feature of a 
neutron star is the value of the maximum mass Mmax for which the star can 
exist for the specific equation of state. As displayed in Fig. 1(c), Mmax is 
strongly dependent on the equation of state, while a stiffer equation leads to 
larger Mmax- 

In Fig. 2(a), we present the information entropy 5", given by Eq. (5), as 
a function of the mass M. We find that S* is a decreasing function of M 
in the region denoting a stable neutron star. The above result is a direct 
consequence of the fact that when the mass of the star increases, its radius 
decreases, so does its volume, while its energy density (or its mass density) 
becomes more localized. Thus, the star is less extended, more compact and 
S is smaller. The effect of the parameter c is just to shift the curve S versus 
M. 

In Fig. 2(b) we plot the information content H — employed in the 
LMC definition C ^ H-D. It is seen that both S{M) and H{M) exhibit the 
same monotonic trend as functions of M. 

In Fig. 2(c) we display the disequilibrium D{M). Increasing M corre- 
sponds to a more concentrated density distribution, its energy density be- 
comes more localized, resulting to a monotonically increasing D. The rate of 
this increase is clearly greater in the region close to the value of Mmax- This 
is due to the fact that as M approaches to its maximum value, it becomes 
almost independent of R. We also observe a reciprocal behavior of the trends 
of S{M) and D{M), as expected from their definitions. 

Complexity C, is plotted in Fig. 2(d). In the region denoting a stable 
neutron star, C is a monotonically decreasing function of the star mass M. 
The most interesting result in this figure is that a neutron star can not grow 
in complexity as its mass increases towards the limit of Mmax- Considering 
the fact that in nature the most probable values for the mass of a neutron 
star vary between 1.4 Mq and 3 M©, we note that a neutron star is eventually 
a physical system of minimum complexity. It is an ordered system, since in 
the corresponding region the rate of decrease of C becomes very small and 
can be considered as a plateau of minimum (zero) complexity. 
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This result becomes more striking in the following set of figures, Fig. 3, 
where we plot in three-dimensions (3D) information and complexity mea- 
sures, as functions of both M and i?, taking advantage of the fact that 
each choice of initial values in the equation of state provides a different pair 
{R, M), reflecting the competition between the gravitational and degenerate 
gas pressures. The facts that the most probable radii of a neutron star are 
close to 10 Km, together with the aforementioned comment on the most 
likely masses, lead us to conclude that a neutron star is in general, a system 
of minimum complexity. Furthermore, it can not grow in complexity as the 
mass or radius increase inside the regions imposed and commented above. 
The neutron star is an ordered system. From the 3D plots of Fig. 3 we can 
visualize the variation of S and D as functions of R, by keeping M constant. 
Information entropy S is an increasing function of i?, i.e. a larger radius 
corresponds to a larger volume, the energy density (or its mass density) be- 
comes more delocalized, and hence the system is more extended and the 
information describing it increases. On the other hand, the disequilibrium D 
is a monotonically decreasing function of R corroborating the fact that the 
system tends to cquiprobability as R increases. 

In order to study in more detail the connection between the information 
and complexity measures with the nuclear interaction and gravity, we plot 

D and C (which correspond to a maximum mass of a neutron star), as 
functions of R and M, varying c and G respectively. In Fig. 4 and Fig. 5 we 
present the results for the effect of the nuclear interaction in two cases. This 
is done by modifying the equation of state by varying c from 0.7 to 1.5 and 
then we plot information and complexity measures, flrst as functions of M 
in Fig. 4 (for a fixed value of R — 11.5 Km) and second, of R (for a fixed 
value of M = 1.5 M©) in Fig. 5. 

Therefore, by keeping R fixed and studying the dependence of S*, D and C 
on the nuclear interaction indirectly employing M, we note that S is almost 
a linear decreasing function of M (Fig. 4(a)), D is increasing exponentially 
(Fig. 4(c)), resulting to a fast exponential decrease of C versus M (Fig. 4(d)). 
Since R is fixed and so does the volume, increasing M corresponds to a more 
localized energy density, hence S decreases with decreasing D. 

The approximate linear and exponential expressions for S{M), D[M) and 
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C(M) are (Fig. 4, = 11.5 Km): 



5 = -7.371M + 2.581, (22) 
D = 7.027 e2-307M ^ 313.550, (23) 

C = 454.949 - 0.001. (24) 

Approximate expressions have been obtained by the application of the 
least squares fitting (LSF) method. We use hnear relations for fitting of the 
form y = C1X+C2, while exponential relations are of the form y = Ci e'^/'^^+Ca. 

On the other hand, we keep M = 1.5 Mq and we see that all measures 
5", D, and C are linearly depended on i?, as c varies. Information entropy S 
increases (Fig. 5(a)), while disequihbrium D decreases (Fig. 5(c)). For fixed 
M an increasing R corresponds to a larger volume. Then, the star becomes 
more extended and accordingly, its energy density becomes less localized. 
Complexity C increases as a result of that delocalization, since the system 
becomes less ordered (Fig. 5(d)). 

The approximate (fitted) hnear expressions for S{R), D{R) and C{R) are 
(Fig. 5, M = 1.5 Mq): 

S = 0.359i? - 12.713, (25) 
D = -96.405i? + 1611.450, (26) 
C = 0.013i? - 0.060. (27) 

We repeat the same series of calculations examining this time the grav- 
itational dependence of information and complexity measures, by varying 
the gravitational parameter G in the range from 0.9 G to 1.1 G, while the 
equation of state is fixed. Our aim is to see how the variation of G affects 
quantitatively ^S" in a neutron star and compare with the results of the previ- 
ous case. These results are presented in Fig. 6 and Fig. 7. A general comment 
is that the two cases are almost equivalent. The trends and the behavior of 

D and C obtained by varying the gravitational parameter G, keeping the 
equation of state fixed, are almost the same with the corresponding trends 
obtained by varying the equation of state for fixed G. 

Specifically, first, for fixed R = 11.5 Km, S decreases linearly with M 
(Fig. 6(a)), D increases exponentially (Fig. 6(c)), resulting to a fast expo- 
nential decrease of C with M (Fig. 6(d)). The energy density of the system 
becomes more localized as G increases for a fixed R (fixed volume). 



11 



The (fitted) approximate linear and exponential expressions for S{M), 
D{M), and C(M) are (Fig. 6, = 11.5 Km): 

S = -7.081M + 2.027, (28) 
D = 362.941 e°-^i2M _ 531 qq?, (29) 
C = 336.821 e-^-^^^^. (30) 

Second, for fixed M = 1.5 M©, S increases linearly with R (Fig. 7(a)), D 
decreases exponentially (Fig. 7(c)), resulting to a linear increase of C with 
R (Fig. 7(d)). The energy density of the system becomes less localized as G 
increases for a fixed M, as a result of the radius increase. 

The approximate linear and exponential expressions for S{R), D(R), and 



C{R) are (Fig. 7, M = 1.5 M©): 

,S = 0.409i? - 13.348, (31) 

D = 41451.805 e-o-^"^-^ - 143.591, (32) 

C = 0.014i? - 0.067. (33) 



Finally in Fig. 8, we present the direct dependence of complexity C on 
the parameters c and G. It is seen from Fig. 8(a) that complexity for a given 
-^max is a decreasing function of the equation of state parameter c (the trend 
is equivalent with the one in Fig. 4(d)), while it increases exponentially with 
the parameter of the gravitational field Fig. 8(b). 

The corresponding (fitted) exponential expressions for C(c), and C{G) 

are: 

C = 0.006 6-^-2°^^ + 0.003, (34) 
C = 5.45 X lO"'^ e^^-^^^^ - 0.001. (35) 



4 Summary 

We present a study of neutron stars from the point of view of information and 
complexity theories. It is shown that the measures of information entropy 

5 and disequilibrium D can serve as indices of structure of a neutron star. 
More specifically, 5" is a decreasing function of the mass of the star, while 
it is an increasing one of its radius. This result is consistent with the fact 
that as a neutron star's mass increases, its radius decreases resulting to more 
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localized energy and mass densities. The disequilibrium D shows an inverse 
behavior. It is an increasing function of the mass and a decreasing one of its 
radius. More localized energy and mass densities correspond to a distribution 
far from equiprobability and as a result the disequilibrium of the system is 
higher i.e. it is far from equilibrium. 

The complexity C of a neutron star is a decreasing function of its mass. 
It almost vanishes for a vast set of pairs of valucs(i?, M), while it increases 
rapidly for masses less than 1.5 Mq and radii greater than 12 Km. But this 
is a not such a favorable case for a neutron star compared with astronomi- 
cal observations done so far. The favorable one, for masses larger than 1.5 
Mq and radii less than 12 Km corresponds to almost vanishing complexity, 
supporting the conclusion that a neutron star is an ordered system, which 
cannot grow in complexity as its mass increases. 

Furthermore, we investigate the impact of the equation of state parameter 
c and the gravitational parameter G on S and C. The behaviors of infor- 
mation and complexity measures are equivalent in both cases. Complexity 
decreases exponentially with the mass, while it increases linearly with the 
radius. In direct calculations, complexity decreases exponentially with the 
equation of state parameter c, while it increases exponentially with the grav- 
itational parameter G. 
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5 Figures 
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Figure 1: (a) Symmetry Energy vs baryon density ra, (b) Energy vs Pressure, 
and (c) Mass vs Radius. 
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Figure 2: (a) Entropy S{M), (b) Information Content H{M), (c) Disequi- 
librium D{M), and (d) Complexity C{M). The insets are commented in the 
text. 
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Figure 3: 3D display of (a) Entropy S{R,M), (b) Disequilibrium D{R,M), 
and (c) Complexity C{R,M), projected for each case on two planes: (a) 
R- M and S - R, (h) R- M and D ~ R, (c) R-M and C - R. 
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Figure 5: (a) S{R), (b) H{R), (c) D{R), and (d) C{R), by varying c for a 
fixed mass M = 1.5 M© (see text, Eqs. (25)-(27)). 
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2.5 




Figure 6: (a) S{M), (b) H{M), (c) D{M), and (d) C(M), by varying G for 
a fixed radius i? = 11.5 Km (see text, Eqs. (28)-(30)). 
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Figure 7: (a) S{R), (b) H{R), (c) D{R), and (d) C{R), by varying G for a 
fixed mass M = 1.5 M© (see text, Eqs. (31)-(33)). 
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Figure 8: (a) Complexity vs the equation of state parameter c, and (b) 
Complexity vs the gravitational parameter G, for a given M = Mmax = 1-5 
M0 (see text, Eqs. (34)- (35)). 
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